Can bimodality exist without phase transition? 
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Abstract 

Here we present an explicit counterexample to the widely spread beliefs about an exclusive 
role of bimodality as the first order phase transition signal. On the basis of an exactly solvable 
statistical model generalizing the statistical multifragmentation model of nuclei we demonstrate 
that the bimodal nuclear fragment size distributions can naturally appear in infinite system 
without a phase transition. It appears at the supercritical temperatures due to the negative 

£T) values of the surface tension coefficient. The developed statistical model corresponds to the 

compressible nuclear liquid with the tricritical endpoint located at one third of the normal 

p^i nuclear density. 

U 

<£ 1 Introduction 

p,i Nowadays the bimodality is considered as an unambiguous signal of the first order phase 

transition (PT) in finite but large systems. The authors of such beliefs [H [21 [3] identify each 
local maximum of the bimodal distribution with a pure phase. For instance, T. Hill justified 
his assumption on bimodality appearance in finite systems by stating that due to the fact that 
an interface between two pure phases 'costs' some additional energy, the probability of their 
coexisting in a finite system is less than for each of pure phases pQ. At the same time it is 

i ", believed [H [21 [3] that in the thermodynamic limit a bimodality corresponds to a mixed phase 

only. 

Here we give an explicit counterexample based on the exact analytical solution of the con- 

f-^ strained statistical multifragmentation model (CSMM) of nuclei in the thermodynamic limit 

which leads to the bimodal fragment distributions inside of the cross-over region without the 

,y~ phase transition existence. In addition, we develop a realistic equation of state for the liquid 



phase which, in contrast to the original SMM formulation [4], is a compressible one. The sug- 
gested approach obeys the L. van Hove axioms of statistical mechanics. The second important 
element of the present model is a realistic parameterization of the surface tension temperature 
dependence which is based on the exact analytical solution of the partition function of surface 
deformations. 

• i-H 

X 

2 CSMM with compressible nuclear liquid in ther- 
modynamic limit 

The general solution of the CSMM partition function formulated in the grand canonical vari- 
ables of volume V, temperature T and baryonic chemical potential fj, is given by [5] 

Z(V, T, M ) = £ e A « V [l - qg^] ~* , (1) 

{An} 

where the set of A n (n = 0, 1, 2, 3, ..) are all the complex roots of the equation 

X n = F(V,\ n ), (2) 



ordered as Re(X n ) > Re{\ n+ \) and Im(Xo) = 0. The function T(V, A) is denned as 

T(uu ( mT \ l \ ti-\Tb \ ^ f ( Pl (T,p)-XT)bk \ 
F(V,\)= I — 1 ziexpj J.+2^ <^ fe (r)exp| J.. (3) 



fc=2 



Here m ~ 940 MeV is a nucleon mass, zi = 4 is an internal partition (the degeneracy factor) 
of nucleons, b = 1/ po is the eigen volume of one nucleon in a vacuum (po — 0.17 fm 3 is the 
normal nuclear density at T = and zero pressure) . The reduced distribution function of the 
fc-nucleon fragment in (I3J) is defined as 



3 



f>k>i(T) = ( — j fc r exp 



a(T) tf 



(4) 



where r ~ 1.825 is the Fisher topological exponent and cr(T) is the T-dependent surface 
tension coefficient. Usually, the constant, parameterizing the dimension of surface in terms of 
the volume is <j = | . 

In ([3]) the exponents exp(— Xbk) (k = 1,2,3,...) appear due to the hard-core repulsion 
between the nuclear fragments [5], while pi(T, p) is the pressure of the liquid phase. 

We consider the thermodynamic limit only, i.e. for V — > 00 it follows K(V) — > 00. Then 
the treatment of the model is essentially simplified, since Eq. (I2j can have only two kinds of 
solutions [5], either the gaseous pole p g (T,p) = TXo(T,p) for ,F(V, Ao — 0) < 00 or the liquid 
essential singularity pi(T, p) = TXo(T, p) for F(V, Ao — 0) — > 00. The mathematical reason why 
only the rightmost solution Xq(T, p) = max{i?e(A n )} of Eq. ((2]) defines the system pressure is 
evident from Eq. (IT]): in the limit V — >• 00 all the solutions of ([2]) other than the rightmost 
one are exponentially suppressed. 

In the thermodynamic limit the model has a PT, when there occurs a change of the right- 
most solution type, i.e. when the gaseous pole is changed by the liquid essential singularity 
or vice versa. The PT line p = p c (T) is a solution of the equation of 'colliding singularities' 
p g (T, p) = pi(T, p), which is just the Gibbs criterion of phase equilibrium. The properties of a 
PT are defined only by the liquid phase pressure pi(T, p) and by the temperature dependence 
of surface tension cr(T). 

In order to avoid the incompressibility of the nuclear liquid we suggest to employ the 
following parameterization of its pressure 

W(T)+p + a 2 (p-po) 2 + a A (p-p )* 
Pi = 1 • ( 5 ) 

Where W(T) = Wo + jy- denotes the usual temperature dependent binding energy per nucleon 
with Wo = 16 MeV, while the constants p = -W , a 2 ^ 1.233 • 10~ 2 MeV^ 1 and a 4 ^ 
4.099 • 10~ 7 MeV~ 3 are fixed by the requirement to reproduce the normal nuclear matter 
properties, i.e. at vanishing temperature T = and normal nuclear density p = po the liquid 
pressure must be zero. 

In addition to the new parameterization of the free energy of the A:-nucleon fragment ([3]) 
we propose to use the following parameterization of the surface tension coefficient 



cr(T) = ao 



T 

± cep 



T 
± cep 



c 

sign(T cep - T) , (6) 



with C = const > 1, T cep = 18 MeV and a = 18 MeV the SMM. In contrast to the Fisher 
droplet model [6] and the usual SMM, the CSMM surface tension ^ is negative above the 
critical temperature T cep . 
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Figure 1: Model phase diagrams are shown in T — /i plane (a)) and p — p plane (b)). In the panel 
a) a first order PT is shown by the solid curves. The vertical dashed lines display the second order 
PTs and the black circles correspond to the tricritical endpoints marked as 1 (nuclear matter) and 
2 (antinuclear matter). A cross-over occurs along the dotted vertical line of the vanishing surface 
tension coefficient. The grey areas in the panel b) show the mixed phases of the first order PTs. 
The isotherms are shown for T = 11, 16, 17, 18 MeV from bottom to top. Negative baryonic charge 
densities correspond to an antimatter. 

Now we would like to study the fragment size distribution in two regions of the phase 
diagram in order to elucidate the role of the negative surface tension coefficient. In order 
to demonstrate the pitfalls of the bimodal concept of Refs. [21 [3j Q] we study the gaseous 
phase and the supercritical temperature region, where there is no PT by construction. As 
one can see from Fig. [2] in the gaseous phase, even at the boundary with the mixed phase, 
the size distribution is a monotonically decreasing function of the number of nucleons A; in a 
fragment. However, for the supercritical temperatures one finds the typical bimodal fragment 
size distribution for a variety of temperatures and chemical potentials as one can see from 
Figf2j 

A sharp peak at low k values reflects a fast increase of the probability density of dimers 
compared to the monomers (nucleons), since the intermediate fragment sizes do not have the 
binding free energy and the surface free energy and, hence, the monomers are significantly 
suppressed in this region of thermodynamic parameters. On the other hand, it is clear that 
the tail of fragment distributions in Fig. [2] decreases due to the dominance of the bulk free 
energy and, hence, the whole structure at intermediate fragment sizes is due to a competition 
between the surface free energy and two other contributions into the fragment free energy, i.e. 
the bulk one and the Fisher one. 



3 Conclusions 



In the present work we showed that the bimodal distributions can naturally appear in infinite 
system without a PT. At the supercritical temperatures a bimodal distribution is generated by 
the negative values of the surface tension coefficient. This result is in line with the previously 
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Figure 2: Fragment size distribution in the gaseous phase is shown for a fixed baryonic chemical 
potential /i = —27.5 MeV and three values of the temperature T. The dotted curve is found exactly 
at the boundary of gaseous and mixed phases. 

discussed role of the competition between the volume and the surface parts of the system free 
energy. 

Also we suggested the new parameterization of the CSMM liquid phase pressure which 
repairs the two main pitfalls of the original SMM and allows one to consider the compress- 
ible nuclear liquid which has the tricritical endpoint at the one third of the normal nuclear 
density. Surprisingly, the suggested approach to account for the nuclear liquid compressibility 
automatically leads to an appearance of an additional state that in many respects resembles 
the physical antinuclear matter. 
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